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- - - Abstract 

(^ ' We consider the leading one-chain term in a skeleton expansion for QCD observables and show that for 

p^ I energies Q^ > A^, where Q^ — A^ is the Landau pole of the coupling, the skeleton expansion result is 

&X), equivalent to the standard Borel integral representation, with ambiguities related to infrared (IR) renor- 

.^ . malons. For Q^ < A^ the skeleton expansion result is equivalent to a previously proposed modified Borel 

fvq ' representation where the ambiguities are connected with ultraviolet (UV) renormalons. We investigate the 
Q -dependence of the perturbative corrections to the Adler-Z? function, the GLS sum rule and the polarised 

^^ , and unpolarised Bjorken sum rules. In all these cases the one-chain result changes sign in the vicinity of 

C^^ . Q^ — A^, and then exhibits freezing behaviour, vanishing at Q^ = 0. Finiteness at Q^ — A^ implies specific 

^S| ' relations between the residues of IR and UV renormalon singularities in the Borel plane. These relations, 

(^ ' only one of which has previously been noted (though it remained unexplained), are shown to follow from 

o 



the continuity of the characteristic function in the skeleton expansion. By considering the compensation 
(^ , of non-perturbative and perturbative ambiguities we are led to a result for the Q^-dependence of these ob- 

servables at all Q^, in which there is a single undetermined non-perturbative parameter, and which involves 



p. 



flj ' the skeleton expansion characteristic function. The observables freeze to zero in the infrared. We briefly 

consider the freezing behaviour of the Minkowskian i?e+e- ratio. 
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I. INTRODUCTION 

Thanks to asymptotic freedom, fixed-order QCD perturbation theory can potentially provide 
accurate approximations to physical observables at suitably large energy scales, Q^ . Such a per- 
turbative description necessarily breaks down below the Landau singularity at Q^ = A^, and the 
infrared behaviour unavoidably involves non-perturbative effects. In fact non-perturbative infor- 
mation is needed even to make sense of perturbation theory, since higher perturbative coefficients 
exhibit factorial growth, and the perturbation series is not convergent. Using a Borel integral to 
represent the resummed perturbation series, the Borel integral is ambiguous due to singularities 
on the integration contour along the positive real semi-axis in the Borel plane, so-called infrared 
(IR) renormalons. These ambiguities are structurally the same as terms in the operator product 
expansion (OPE) in powers of h? /Q'^. OPE ambiguities and Borel representation ambiguities can 
compensate each other, allowing the perturbative Borel and non-perturbative OPE components to 
be separately well-defined once a regulation of the Borel integral, such as principal value (PV), has 
been chosen [jJl. For Q^ < A"^, however, the Borel representation which is correlated with terms 
in the OPE breaks down. In a recent paper Ref.[2], which focussed on the infrared freezing of the 
Minkowskian R^+e- ratio, it was suggested that below Q^ = A^ one should use a modified Borel 
representation whose ambiguities come from singularities lying on the integration contour along 
the negative real semi-axis, so-called ultraviolet (UV) renormalons. This Borel representation has 
ambiguities which are structurally the same as a modified expansion in powers of Q^/A^, and once 
regulated both components can remain defined in the infrared. This change of Borel representa- 
tion has been claimed not to be physically motivated in Ref.p|, where different conclusions about 
infrared behaviour are reached. In this paper we shall show that if we postulate a QCD skeleton 
expansion [J, |5| , then the leading one-chain term reproduces the standard Borel representation for 
q2 > ^2^ ^^^ ^]^g proposed modified Borel representation for Q^ < A'^. 

We consider the infrared behaviour of the one-chain result for some Euclidean QCD observables. 
We shall concentrate on the Adler-D function, the GLS sum rule and the polarised and unpolarised 
Bjorken sum rules jaild]- The skeleton expansion result automatically freezes to zero as Q^ -^ 0. 
For the observables we consider, the freezing to zero occurs after the Borel resummed perturbative 
corrections to the parton model result change sign in the vicinity of Q^ = A^. Individual renormalon 
contributions to the Borel integral diverge at Q^ = A^ , but we find that when all of the renormalons 
are summed over, one obtains a finite result. This finiteness requires relations between the residues 
of infrared and ultraviolet renormalons. Only one of these relations has previously been noted 
[^ , and we show that they arise from the continuity of the characteristic function in the skeleton 
expansion. Considering the compensation of perturbative and OPE ambiguities alluded to above, 
we are led to an expression for the Q^-dependence of the observable written in terms of the 
characteristic function, and containing a single undetermined non-perturbative parameter. This 



result freezes to zero in the infrared. Existing discussions of infrared freezing behaviour have largely 
focused on the Analytic Perturbation Theory (APT) approach 9]. In this formalism one expands 
observables in a basis of functions which have smooth infrared behaviour. For Euclidean observables 
the unphysical Landau singularity in the coupling is cancelled by a power-like correction. In 
contrast in our discussion finiteness and continuity emerge thanks to a subtle interplay between 
UV and IR renormalons. 

The plan of the paper is as follows. In Section 2 we shall introduce the QCD skeleton expansion, 
and show that the one-chain leading term is equivalent to the standard Borel representation for 
q2 y ^2^ ^^^ |.Q ^YiQ modified representation for Q^ < A'^. We discuss what can be learnt about 
the infrared freezing of observables. In Section 3 we describe the Borel plane renormalon structure 
for our chosen Euclidean observables, and we show that finiteness at Q^ = A^ only holds if there 
are cancellations between the residues of IR and UV renormalons, the cancellations rely on a 
previously unknown relation between the IR and UV residues. We write down a result for the 
Q^ dependence of the resummed observables in terms of Exponential Integral (Ei) functions, and 
plot the infrared freezing behaviour to zero noted above. In Section 4 we consider the skeleton 
expansion for the Adler D function, and give an expression for the characteristic function of the 
leading one-chain term. Making a power series expansion, and changing variables, we explicitly 
obtain the Borel representations, and relate the IR and UV renormalon residues to the power 
series coefficients of the characteristic function. Continuity of the characteristic function is shown 
to underwrite the relations between UV and IR renormalon residues noted above. In Section 5 
we derive the result for Q^-dependence including non-perturbative effects mentioned above. In 
Section 6 we briefly consider Minkowskian observables, specifically i^g+g-, and modify some of the 

n n 

conclusions of Ref.';2!] in the light of the criticisms of Ref.|3|. Section 7 contains a discussion and 
our conclusions. 



II. QCD SKELETON EXPANSION AND BOREL REPRESENTATIONS 

Consider a generic Euclidean QCD observable ViC^) having the perturbative expansion 

Vpt{Q^) = a{Q^) + Y,dna^^\Q^) . (1) 

n>0 

Here a{Q'^) = as{Q'^)/ir is the renormalised coupling. Throughout this paper we will use the 
one-loop approximation for the coupling, 

where b = (33 — 2Nf)/6 is the leading beta- function coefficient in SU(3) QCD with Nf active 
quark flavours. Q^ = —q^ > is the single spacelike energy scale. As Q^ ^ oo asymptotic freedom 
ensures that T){Q^) — > 0. Our interest is in the infrared limit Q^ — > 0, and the infrared behaviour 
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FIG. 1: A chain of fcrmion bubbles with momentum k running through them. 



of 'D{Q'^). Specifically, is it possible that freezing to a finite infrared limit 1^(0) occurs? This is 
an intrinsically non-perturbative question which cannot be answered by perturbation theory alone. 
One has in addition the non-perturbative contribution arising from the operator product expansion 
(OPE), 



Vnp{Q^ 



^^\q^ 



(3) 



The freezing limit, if any, of T>{Q'^) = 'DpxiQ^) + T^np{Q'^), depends on the behaviour of both 
components as Q^ ^ 0. Perturbative freezing will not arise from fixed-order perturbation theory, 
one needs an all-orders resummation of Eq.l^Q). Unfortunately our exact information about the 
higher-order coefficients is limited, at best, to calculations of di and d2, higher-orders are unknown. 
All-orders information is only available in the large- A^j limit where one expands each d„ as 



dn 



dt^m + dt~^^m 
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/ 



+ ...+ 



dM. 



(4) 



The leading large- A^/ coefficient dn can be computed exactly to all-orders since it arises from a 
restricted set of Feynman diagrams in which a chain of n fermion bubbles (shown in Fig. ^ is 
inserted in a basic skeleton diagram [lO|, lllll. In principle one can consider more than one chain 
and construct a QED skeleton expansion [l2]. In QCD one can replace A^^ by (33/2 — 36) , and 
obtain an expansion in powers of 6, 



d„ = d(«)6" + 4"-i)5"-i + . . . + 4°) 



,{L) 



An) 



0.0 



(5) 



and an all-orders 



The leading-5 term dn ' = dn b"' can then be used to approximate dn 

resummation of these terms performed to obtain 'Dpj,{Q'^). Use of the one- loop form of the coupling 

in Eq.@ ensures that this resummed result is RS-independent. 

The leading term of the skeleton expansion arises from integrating over the momentum k flowing 
through the chain of bubbles [J, [sL ll^ 



v^pkQ' 



dtLo{t)a{e~'tQ 



(6) 



Here t = k"^ /Q"^, and uj{t) is the so-called characteristic function of the observable. The constant C 



depends on the subtraction procedure used to renormalise the bubble. Standard MS subtraction 
corresponds to C = —5/3. ^From now on we shall assume C = which corresponds to the so-called 



^-scheme, MS subtraction with renormalization scale fi^ = e ^"Q^. A in Eq.@ will refer to that 
in the ^/-scheme. The characteristic function satisfies the normalization condition 



dt u;{t) = 1 , (7) 

'0 

which ensures the leading a{Q'^) coefficient of unity assumed in Eq.(^. The form of Lo{t) changes 
at t = 1, and the range of integration splits into an IR and a UV part 



T^pW) = dt uiR{t)a{tQ'^) + dt uJuv{t)aitQ^) , 



(8) 



the IR part corresponding to k"^ < Q^, and the UV part to k'^ > Q^. By making a change of 
variable one can transform the leading skeleton term into a Borel representation. For Q^ > A^ 
one has the standard Borel representation (we shall explicitly write down the required changes of 
variable in Sec. 4), 

Jo 
Here B['D]{z) is the Borel transform, defined by, 

B[V^,%z) = Y^'-^. (10) 

n=0 

B[T>pj]{z) contains singularities along the real z-axis. In the large-6 approximation these are 
single and double poles at positions z = z„ and z = —Zn, with z„ = 2n/b, n = 1,2,3.... The 
singularities on the positive real semi-axis are referred to as infrared renormalons, IRn, and those 
on the negative real semi-axis as ultraviolet renormalons, UVn- The IRn renormalons cause the 
Borel representation to be ambiguous since they lie on the integration contour along the positive 
real z-axis. The difference between routing the contour above or below the singularity yields an 
ambiguity 



A^S~{^) > (11) 



A2X" 

which has the same form as a term in the OPE in Eq. (jSJ , so that OPE ambiguities associated with 
the (A^/Q^) OPE term in T>f^p{Q'^) can potentially cancel against the IRn renormalon ambiguity 
allowing each component separately to be well defined [l|. In practice we shall choose to take a 
Principal Value (PV) definition of the integral. The IR part of the t integration in Eq.® produces 
the IR renormalon part of the Borel representation, and needs to be PV regulated. The second 
UV component produces the UV renormalons and does not require regulation. As we shall see in 
the next section the standard Borel representation of Eq.® for Euclidean quantities diverges like 
lna((5^) at Q^ = A^ for each individual IRn or UVn renormalon contribution. When the full set is 
resummed, however, the In a divergence is cancelled and a finite result is found. We shall explore 
this further in Sections 3 and 4. 



For Q^ < Is? , a{Q?) < 0, and the representation of Eq.® is invalid. The key point is that the 
change of variable from t to z is proportional to a((5^), and so if a{Q'^) changes sign the limits of 
integration in z change sign, yielding the modified Borel representation 

/•— oo 

vf^^Q^)= dze-'/''^'^'^B[vfi]{z). (12) 

JO 

This is the modified Borel representation proposed in Ref.^] where it was motivated as a stan- 
dard Borel representation corresponding to an expansion in |a((5^)| = —a{Q'^), since by changing 
variables one can write Ea. (|12|) as 

/■oo 

V^jt^iQ^) = - / dze-^/l'^(«')ls[pi,''^](-z) . (13) 

Jo 

So we see that the one-chain skeleton contribution of Eq.@ is equivalent to the standard Borel 

representation of Eq.Q for Q^ > A^, and to the modified representation of Eo. 1)12(1 for Q^ < A^. 

Note that when we substitute Eq. ()1U() into the Borel representation of Eq. (|13|) . then it reproduces 

the correct form of the perturbative expansion in Eq. (^, for negative a. The modified Borel 

representation now has a contour of integration along the negative real semi-axis, and so it is 

rendered ambiguous by the ultraviolet UVn renormalon singularities. Correspondingly the IR 

component of Eq.Q is now well-defined and it is now the UV component which requires regulation. 

The ambiguity from routing the contour is now 

APi§(Q2)~(|J)". (14) 

It was suggested in Ref.jj] that the usual OPE of Eq.® breaks down for Q^ < A^, as does 
the associated PT Borel representation of Eq.®, and should be recast and replaced by a modified 
expansion in powers of Q^/A^, 

'DNp{Q')=J2Cn(%Y ■ (15) 

The nth term in this expansion has then structurally the same form as the ambiguity associated with 
the UVn renormalon contribution. It was further suggested in Ref.t2i that a Cq term independent 
of Q^ could arise from rearrangement of the standard OPE. This was motivated by a simple 
toy example. In fact in the one-chain approximation no such term arises and both PT and NP 
components freeze to zero. The terms in Ea. (|15() are then in one-to-one correspondence with the 
UVn renormalon ambiguities. From its definition, the QCD skeleton expansion implies Ppt(O) = 
in the Q^ ^ limit. For the one-chain term in Eq.® this simply follows because as Q^ — > the 
integrand vanishes everywhere in the range of integration, since a{tQ?) -^ for any given t. Higher 
multiple chain terms will contain products of the form a{tiQ'^)a{t2Q'^) ... in the integrand and will 
similarly vanish. This then implies that in the infrared limit 'Dnp{Q'^) behaves as 

VMP{Q^)^k(^) , (16) 






FIG. 2: Leading large-Nf contributions to the vacuum polarisation function at nth order in perturbation 
theory. 

where UVn^ is the UV renormalon singularity nearest to the origin in the Borel plane. We should 
note that the modified Borel representation, its infrared behaviour and its connection with UV 
renormalons, has also been discussed in Ref.[la|- In Appendix B of that paper the infrared freezing 
of the Adler function D{Q^) was discussed and it was concluded that from general arguments of 
non-perturbative spontaneous chiral symmetry breaking in the limit of a large number of colours, 
Nc, one expected that as Q^ -^ 0, D{Q'^) -^ like 



D{Q'') 



M2 ' 



(17) 



where M is the mass of a one-meson state, these states remaining massive in the chiral limit. A 



similar result is obtained in Ref.|17t|. Since UVi is the singularity nearest the origin for the Adler 
function, no = 1, and the freezing expectation is indeed consistent with Ea.H16|). Notice that 
strictly the leading behaviour as Q^ ^ is the logarithmic freezing to zero of a{Q'^) contributed by 
the PT component. It is the non-perturbative effects which reflect the UV renormalon structure. 



III. Q2.j)epe]vjj3e]sjce of the EUCLIDEAN OBSERVABLES 



We begin by defining the three Euclidean observables we shall consider. The QCD vacuum 
polarization function, Il{Q'^), is the correlator of two vector currents in the Euclidean region. 



{q^.q. - 9,.uq^MQ^) = Ifivr^i / d^x e^'?-^(0|r[j^(x)>(0)]|0) 



(18) 



The leading-A^j component of Il{Q^) can be calculated from the diagrams in Fig. [3 The Adler 



function, D{Q'^), is then defined via the logarithmic derivative of Ii{Q'^) 

D{Q') = -lQ'^2^iQ')- (19) 

This can be spht into the parton model result and QCD corrections, 2?(Q^), 

(20) 



D{Q') = AT, ^ Q2 / 1 + ^CfV{Q^)^ 



f 

where Nc is the number of colours, Cp = 2N ' ^^^ Qf ^^ ^^^ charge of quark flavour /. Here 
T^iQ"^) = ^pt(Q^) + V iyp JQ'^), with the two components defined as in Eqs. (0) and ©• The 
polarised Bjorken (pBj) 18 1 and GLS [l9| sum rules are defined as 

KpBj = / gf~''''^{x,Q'^)dx 
Jo 



9A 



9V 



1-^CfIC{Q^)\ , (21) 



t> Jo 



'0 

3 



= \l-^CFlCiQ')\ . (22) 

/C((5^) being the QCD corrections to the parton model result, again split into PT and NP compo- 
nents as for D((5^). We have neglected contributions due to "light-by-light" diagrams - which when 
omitted render the perturbative corrections to Kgls and KpSj identical. Finally, the unpolarised 
Bjorken sum rule (uBj) [20| is defined as 



UuBj ^ [' F^''-'nx,Q^)dx 
Jo 



= Ii-^CfU{Q')] . (23) 

The QCD corrections to the parton model result are again split into PT and NP components. 
The leading-A'^j contributions to these three sum rules can be calculated from the diagrams in 
Fig. 131 These large- A^/ results can be used to compute leading-6 all-orders resummations for these 
observables, Vpj.{Q'^), }Cprf,{Q'^) and Up^{Q'^), as described in Section 2. 

■)Wu \ _ v^ Ao{n) - Ai{n)zn , Ai{n)z„ 



The Borel transform of 'Dpl{Q'^) is well-known and can be found in Ref. 



iJ[pig](.) = ^ -°7— \7-- + 



-1 11 + ^) (i + t 



Bo{n) + Bi{n)zn Bi{n)z„ 



n=l (l ^^ 



Zr. 



FIG. 3: Leading large- A^/ contributions to the DIS sum rules of Eqs. (|^ - f^ at nth order in perturbation 
theory. 



Here 



An(n) 



3 n2(n + l)2(n + 2)2 ' 

Boil) = 0, Bo{2) = 1, 
Si(l) = 0, Bi(2) = 



A^(n 



Bo{n) 



4, Biin) 



8 6(-l)"+Hn + |) 
3n2(n + l)2(n + 2)2 

^o(~^) ?i- > 3 
-^i(-n) n>3 



(25) 



. The purpose of the sHght change 



These definitions coincide with Ref.|a|, except for Bi{2) — ^ 

of definition is to make more explicit the single and double pole structure. The Borel transforms 

of /Cpj.('3^) ^^'^ ^pT^Q"^) '^^'^ ^^ found in Refs.ja, I3I, respectively. They have a much simpler 

structure than that of the Adler-D function since they arise from insertion of the chain of bubbles 

into a tree-level diagram, rather than into a quark loop, as shown in Fig. 13 There are only a finite 

number of single poles and no double poles. Consequently we can write out their Borel transforms 

explicitly 

4/9 1/18 8/9 5/18 



S[4^^](z) 



1 + 



and 



Wi 



B[U'pi]{z) 



1/6 



1 + 



+ 



+ 



4/3 



1 



1 



(26) 



1/2 



1 + 



22 



1 



21 



1 



(27) 



As noted in Refs.p, 1^ the leading-6 approximations for the NLO and NNLO coefficients for these 
observables are in reasonable agreement with the known exact coefficients. 

We can now evaluate the Borel integral of Eq.@ to obtain Vpj.{Q'^), lCpj,{Q'^) and Upj,{Q'^). 
Using the integrals 

Jo (1 + z/z„ 



-Zne'-/''Eii-zJa) 



f 

Jo 



dz 



Q-z/a 
{1+z/ZnY 



1 + ^e-^"/'^Ei(z„/c 



(28) 
(29) 



the following resummed expressions are obtained, 

OO f 



n=l 



aiQ' 



a{Q'' 



{AQ{n) - ziAi{n)) - z„^i(n) 



+ {AQ{n) - ZnAi{n)) 






Bo{n) + ZnBiin)) \ , 



^ {Bo{n) + ziBiin)) - z„5i(n) 



(30) 



/cS(Q^) = ^ 



«^^^^"^^^^^<-J^)+-^^^"^^^^^<-^) 



+ lgg-^i/a(Q^)Ei 



^1 ] _ lOe-^^/°(Q')Ei ( -^^ 
a{Q^)) \a{Q') 



(31) 



4?(Q^) = I 



gg-^i/a(Q2)g. 



^1 



a(Q2)y U(Q') 



2e-M«^)Ei(--5^ 
a(Q^J 



Where Ei(x) is the exponential integral function defined (for x < 0) as 

e" 



Ei(x) = 



dt- 



t 



(32) 



(33) 



and for x > by taking the PV of the integral. It has the expansion 

E-i{x) = ln\x\ + -fE + 0{x) , (34) 

for small x, where je = 0.57721 ... is the Euler constant. 

A crucial point is that the above expressions for the Q^-dependence apply at all values of Q^. 
For Q^ < A^ the modified Borel representation, written as an ordinary Borel representation for an 
expansion in powers of \a\, as in Eq. (fT3|) . corresponds to changing a{Q'^) — > —a{Q'^) , z„ -^ —Zn, and 
adding an overall minus sign in Eas. (|30|) - ()32() . One can easily see that these equations are invariant 
under these changes. In Ea. ()30() one needs to change Ai — > —Ai and Bi -^ —Bi, since they contain 
a hidden z„ factor in their definitions, also in Eqs. ipTTj) and (jS21), the prefactor proportional to 1/6 
also needs to change sign since it has been factorised from zi,Z2- The Ei(z„/a(Q^)) functions 
exhibit a logarithmic divergence as their argument goes to zero, and so it would appear that one 



10 



does not obtain a finite result at Q^ = h?. Using Ea. (|5H) . one has, 

2n 



Ei 



6a(Q2 



Ei 



~ Ei 



2 /a2\ 



nlog(QVA 



4%-^ 



IE + ln 



" 1^-' 



(35) 



for A^ ~ Q^. Note that the only terms in Eas. H3U|) - ()32() which could possibly contribute to the 
divergence are e^^"'"Ei(=pz„/a) terms and, as can be seen from Eas. (|28|) and H29() . these are 
generated exclusively by the single pole terms in the Borel transform. The double pole terms only 
generate finite contributions at Q^ = A^. 

Using Eq.jSSl) we obtain the Q^ ^ a2 lij^it of V^p^{Q'^) 



T^mQ") = -Y.'l^Mn) + B,{n)]\n 



n=l 

oo 



»S- 



(36) 



Q' 



+ Y^izlil + JE){-Ai{n) - Bi{n)) + z„(Ao(n) - Bo{n))] + O I ^ - 1 

n=l \ 

So the coefficient of the divergent log term in 'Dp^{Q'^) is, 

oo 
n=l 



(37) 



and for ]Cpj.{Q'^) and Uprf {Q"^) the equivalent coefficients are (—8 + 2 + 16 — 10 = 0) and (8 — 6 — 2 = 
0), respectively. Cancellation clearly occurs in the cases of /Cp^(Q^) and Uprp{Q^) and in the case 
of T>pj.{Q'^) the previously unnoticed relation 



ensures that 'Dpji.{t^) is finite 



A similar relation 



zl+2.Biin + 3) = -zlAi{n) , 



Y,zl[Ai{n) + B,{n)]=d. 

n=l 



Ao(n) = -5o(n + 2) 



(38) 



(39) 



(40) 



was noted in ^8]. We shall show in the next section that the relations of Eas. H38|) and (|40() are 
underwritten by the continuity of the skeleton expansion characteristic function un {t) and its first 
derivative at t = 1. The form of the perturbative corrections, Vp^{Q'^), lCpj.{Q'^) and Up^{Q'^), 
are shown in Fig. [l] 
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FIG. 4: Q^-dependence of the perturbative corrections to the observables in Eas. l|2U|) - (|23|l . rcsummed to all 
orders in the leading-6 approximation. 

Although we have shown that when summed to infinity Eq. (|3U|) is finite at Q^ = A^ , we obviously 
can only plot the expression including a finite number of terms in the n sum. The expression can 
remain finite, however, if we sum the UV renormalons to finite n = N and the IR renormalons 
to n = A^ + 3. In this case the relation of Ea. (|38() will ensure that the divergent terms cancel. 
We took A = 50 and assumed Nf = quark flavours, avoiding the need to match at quark 
flavour thresholds, since we are only interested here in the form of the freezing behaviour, not in a 
phenomenological analysis. 

The plots in Fig. |3 demonstrate two important points about the Euclidean quantities we are 
considering. Firstly the finite behaviour at Q^ = h?, and secondly that the Borel resummed 
perturbative corrections to the parton model result change sign just below or above this point. 
For D these corrections become negative but crucially the full observable D{Q^) remains positive 
at all values of Q^. They then freeze to zero as noted in Sec. 2. 
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FIG. 5: Light-by-light scattering diagrams, used to calculate cou- 
The relation of Ea.H38|) simplifies the expression for the finite part of Ea. ()36() . it becomes 

oo oo 



m 



n=l 

0.123625 . 



n=l 



(41) 



The values /Cpj.((5^ = A^) and Uprf,{Q'^ = A^) are given by a formula identical to Ea.(|IT|). but 
using values of Ao,i(n) and i?o,l(?^) appropriate to /C and U. Although we have not given these 
values explicitly, they are of a much simpler form than in the case of P, and they can easily be 
deduced by comparing Eas. (|26j) and (|27|) with Ea. (|24l) . From this we obtain. 



/cSi^^(Q2 = A2) = ___ln2 



8 
96 



Z^},^^(Q2=A2 



3b 



ln2 . 



(42) 



IV. SKELETON EXPANSION AND BOREL REPRESENTATIONS FOR THE ABLER 
FUNCTION 

We begin with the one-chain skeleton expansion result for the vacuum polarization function 
n(Q2) defined in Eq.dlHl), 



n(g2)= / dt uJuit)aitQ^) , 
Jo 

where the characteristic function wn(i) is given by 

tE{t) t<l ^ IR 



(43) 



a;n(*) 



4 

- < 
3 



(44) 






t>l ^ uv 



It can be obtained from the classic QED work of Ref.j2ll| by simply including appropriate colour 
factors.^ In this language it is related to the Bethe-Salpeter kernel for the scattering of light-by- 
light, and is the first term in a well-defined QED skeleton expansion 12]. The diagrams relevant to 



^ The origin of the minus sign in Eg. 14411 is the difference between the definitions of 11 given in Eq. HSU and Ref.|21 
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offy 



the kernel are shown in Fig. [51 It is easy to see how, by connecting the ends of the fermion bubble 
chain in Fig. ^to the niomentuni k external propagators in Fig. |3 one can reproduce the topo^ 
of the diagrams in Fig. |2 The existence of the QCD skeleton expansion is more problematic 
H(t) is given by |2 



1^. 



(i) = -|l-lnt+(---lnt)t+^ ' 



.2 2 
where L2{x) is the dilogarithmic function. 

L.ix) = 



[L2(-t)+lntln(l + t)]} , 



- In(l-y) 

dy 

y 



(45) 



(46) 



Though we define u;n(i) separately in the IR and UV domains, the two regions are related by the 
conformal symmetry t 



t ■ 



The Adler-D function, related to n(Q^) through Ea. H19|l . will have the one-chain skeleton 
expansion term with characteristic function ujx>{t), 



^S(Q') 



dt ux>{t)a{tQ'^ 



(47) 



coi)it) is obtained from uJu{t) by performing the differentiation of Ea. H19|) on Ea. H43|) and then 
performing integration by parts on the resulting expression. 



^ dQ- JO 



^'piiQn = -iQ' 



dtuju{t)t\ 



_3_ 2^ /■°° d_ 
+ 26^ W Jo dt 



UJuit)t 



3 r 

'^Jo 



dt 



(^u{t) +t—UJu{t) 



ln[a(tQ2)] 
aitQ^) . 



(48) 



The transformation from 11 to 2? therefore induces a transformation in uJu{t) of 



n(Q2 



^W'^^^'^ 



MQ') 



d 4 

uju{t) -^ u;u{t) + t—ujn{t) = --ujv{t) 
dt 3 



(49) 



This transformation spoils the conformal symmetry present in LOu{t)- Indeed the expressions for 
ujx>{t) in the UV and IR regions are slightly more complicated. 



u^l^it) 



-\ I - - In t j t + (1 + t) I L2{-t) + In t ln(l + t) 



^ 1 + lnt+^ + llntU 



(50) 



+ (1 + t) L2(-t"')-lniln(l + t-^) 



(51) 



14 



However, a partial symmetry remains in ujx){t) and this will be elucidated upon in the following 
discussion. We shall now convert the skeleton expansion form into the Borel representations of 
Eqs.@ and ()12|) by making a change of variables. To achieve this it is necessary to write touit) as 
an expansion in powers of t. This yields expressions in both the IR and UV regions comprising an 
expansion plus an expansion times a logarithm. 



.IR/ 



^uit) = - 3 E ^"*" + 1" * E ^""*" • (52) 

\?i=l n=2 / 

The conformal symmetry expressed in Eq. H44|l means that the UV part can also be written in terms 
of the coefficients ^n and ^^ 

, / oo oo \ 



^rW = -:^{^Ut-''-lntJ2Lt- 



(53) 



\n=l 



n=2 



From Eq. (|1S|) , ^„ and ^„ are found to be 

4 (2-6n2)(-l) 



Cn 



>i 



3(n-l)2n2(n + l)2' 



e 



n>l 



2(-ir 



3 (n — l)n{n + 1) 

6 = 



6 = 1 

Performing the transformation in Ea. (|49() allows us to write ujvit) as a similar expansion 

oo oo 



n=l 

oo 



n=2 

oo 



^rw = E[^"(i-^)-^~"]*""+i^*E^'"(^-i)*"" 



(54) 

(55) 
(56) 



n=l 



n=2 



(L)tn2\ 



Using the expansions of Eqs. (|55j) and (|56j) we can now represent 'Dpj.{Q'^) in terms of a Borel 
integral. We take 'Dpj.{Q'^) expressed in terms of ujx>{t) and then split the integral into IR and UV 



regions 



pg(Q2) = / dtiUj,{t)a{tQ^) 
Jo 



oo „i 

E«(Q')y^ dtu'j'{t)\ 



k=0 

oo 



' — ' poo 

+ E"(^')/ rfi^r(t)i 



ha{Q^) 



Int 



Int 



KQ^)E 



5a(Q2) 



fc=0 



„l / OD oo \ 

/ dt V[?n(l + n)+en](t)"+lntVen(n + l)(t)M(lnt)^ 

■^0 \n=l n=2 J 

+ / dti Y^[^^(l-n)-L]{tr-+lntY,Un-mr']ilnt)'' 



,n=l 



n=2 



(57) 
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Where we have used 

"" _ ba{y) 

k=0 



a{xy) = a{y)J2[ "^^Inx . (58) 



We note that [^n(l ~ n) — (,n] =0 for n = 1, which ahows us to omit this term from the above 
sum. This expression may be transformed into a Borel integral of the form of Eq.® by changes 
of variables and integration by parts. We use the change of variables z = —a(Q^){n + 1) Int and 
z = a{Q^){n — l)lnt for IR and UV parts, respectively. Integration by parts is necessary for the 
integrals with an extra Int term. For Q^ > A^, a{Q'^) > 0, we then obtain the standard Borel 
representation, of Eq.® 

•^0 Ln=l "^^ '-- 2(^1+1} n=2^^^'-) (1 



2 ^ b. ^' 

2(n+l) 



+/ dz e"(«') 
Jo 



oo 



[U^-n)-^n] 1 , ^Un-1) 1 



n-1 1+ ^bz Z^ (n-l)2 /^ , 5^ n2 

n=2 + W=T) „=2 ^ W 1 + 2(^ 



(59) 



and for Q^ < A^, a{Q'^) < 0, we obtain the modified Borel representation of Ea. (|12() . in which 
the upper limit in z is — cxd. Having obtained the Borel transform we can now make contact with 
Ea. ()24() and this allows us to make the identifications 

^^^^If^^^ = -Bi{n + l)zn+i n > 1 (60) 

gn(l -n)-in 

n — 1 
for the single pole residues and 

L{n + 1) 



(n + l)2 
Lin - 1) 



Ai{n-l)zn-i n>2 (61) 

Boin + l) + Bi{n + l)zn+i n>2 (62) 

Ao{n-l)-Ai{n-l)zn^i n>2 (63) 



(n-l)2 

for the double pole residues. Substituting the form of ^„ and ^n given by Eq. (|54j) . and comparison 
with Eq. (|25j) . verifies the above equations. 

Equations. (|^n|) - (|H5|) can be used to rewrite the uj^{t) and iJj^ expansions of Eas.lj^ and 
(|56|) in terms of the ^o(^)) ^i(^)) ^-nd BQ{n), Bi{n) renormalon residues. One finds 

, oo oo 

^v%) = 2 E -^l+iBi{n + l)e - Int ^ (n + iflB^in + 1) + z„+ii?i(n + l)]t" (64) 

n=l n=2 

, oo oo 

^v'it) = 2 E ^n-i^i(" - 1)*"" + Int ^ (n - l)'[Ao(n - 1) - z„_iAi(n - 1)]*"". (65) 

71=1 n=2 
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The discontinuity at t = 1 is then found to be 



oo 



^g^(l) - u:i,^{l) = I ^ zl[A,{n) + B^{n)] , (66) 



2„=i 



which vanishes using Ea. iP^ . In the language of ^„ and ^„ coefficients, Ea.(|HH]) is equivalent to 

oo 

-2^ (<„ + !„) =0. (67) 



ra=l 



So the relation between UV and IR renormalon residues of Eq.(j2Hl), which guarantees finiteness at 
Q^ = A^, ensures that the characteristic function uj-pit) is continuous at t = 1. 
For the first derivative at t = 1 one finds the discontinuity 



dw^^ 



dt 



do,^^ 



dt 



t=l n=l 

,2 oo ,2 oo 

- -Y^zl[Ao{n) + Bo{n)] + -Y,4[Mn)-B,{n)], (68) 

n=l 71=1 

Equation. H39|) . which ensures that I?p^(A^) is finite, means that the first line of this expression 
vanishes. The second line also vanishes, ensuring continuity of the first derivative of ujx>{t)- This 
also ensures that the 'Dpj. (A^) is finite (the prime denoting the first derivative d/dlnQ). Indeed, 
the required relation corresponding to the vanishing of the coefficient of the potentially divergent 
In a term in 'Dpj. (A^) is, 

oo 

Ypzl{M{n) - Bi{n)) - zliAoin) + Bo{n)] = . (69) 

n=l 

So finiteness of the first derivative of V^ '{A) at Q = A, corresponds to continuity of the first 
derivative of uj{t) at t = 1. Furthermore, Ea.H69|) written in terms of ^n and ^„ is simply Ea. H67() . 
with an extra factor of —2. Consequently, the continuity of ujx>{t) and its first derivative stem for 
a single relation, Eq. 1)6 7|) . The second and third derivatives are also continuous at i = 1, and their 
discontinuities involve additional new structures built from combinations of the ^o,i and i?o,i- To 
ensure finiteness of l^pj' (Q^) at Q^ = A^, one requires the relation 

oo 

^[34(^i(n) + Bi{n)) - 24(Ao(n) - So(n))] = . (70) 

n=l 

For finiteness of V^^' (Q^) at Q^ = A^ one requires the relation 

oo 

5][44(^i(n) - Bi{n)) - 34(^o(n) + Bo{n))] = . (71) 

n=l 

Eas. lfTn]) and (f7T|) are also required in order for the second and third derivatives of oJvif) to be 
continuous at t = 1, furthermore, they can both be derived from the following relation 

oo 

^ (n3e„ + Sn^^n) = . (72) 

n=l __ 

The fourth and higher derivatives of oj-d (t) are discontinuous at t = 1 as noted in [J] . 
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V. SKELETON EXPANSION AND THE NP COMPONENT 

In this section we wish to consider more carefuhy the compensation of ambiguities between 
renormalons and the OPE. The regular OPE is a sum over the contributions of condensates with 
different mass dimensions. In the case of the Adler function the dimension four gluon condensate 
is the leading contribution, 

Go(a(Q2)) = ^iS)\GG\Q)CGG{a{Q^)) , (73) 

where Cgg{(^{Q'^)) is the Wilson coefficient. In general the n term in the OPE expansion of 
Eq.Q will have the coefficient 

C„(a(Q2)) = Cn[a{Q^)f"(.l + 0{a)) . (74) 

The exponent 5n corresponding to the anomalous dimension of the condensate operator concerned. 
Non-logarithmic UV divergences [24] lead to an ambiguous imaginary part in the coefficient so that 
Cn = Cn ± iCn ■ If One considers an IRn renormalon singularity in the Borel plane to be of the 
form Kn/{1 — z/zn)'^" then one finds an ambiguous imaginary part arising of the form 

lm[VpT] = ±K„^^e-^"/'^(Q')ai-T" [1 + 0{a)] . (75) 

r(7n) 

Here the it ambiguity comes from routing the contour above or below the real z-axis in the 
Borel plane. This is structurally the same as the ambiguous OPE term in Eq. (|74|) . and if Cn = 
KniTz2" /T(yr,) and 5n = 1 — jn, then the PT Borel and NP OPE ambiguities can cancel against 
each other 23] . Taking a PV of the Borel integral corresponds to averaging over the it possibilities. 



For Q^ < A^ the modified expansion of Eq. ()15() will have an n**^ coefficient of the form 

Cn{a{Q^)) = Cn[a{Q^)f{l + 0{a)) . (76) 

Now the exponent 5n is related to the anomalous dimension of dimension 6, four-fermion op- 
erators associated with UV renormalons [2J|, IR divergences associated with these render the 
imaginary part ambiguous, and Cn = Cn ± Cn ■ The modified Borel representation of Ea. H12|) 
has ambiguities arising from UV renormalons. Assuming that the UVn singularity is of the form 

kn/{l + z/znf" one finds 

Im[PpT] = ±^„— ^e^-Z'^^Q )ai-^" [I + 0{a)] . (77) 

i(7nj 

This is structurally the same as the ambiguity in the modified NP expansion coefficient in Ea. H76() . 
and if Cn = KnT^z]^ /T{'^n) and ^n = 1 — 7n, the ambiguities can be cancelled. 

In the one-chain (leading-6) approximation the renormalons are single or double poles corre- 
sponding to 7 = 1 or 7 = 2, and correspondingly the ambiguous imaginary parts in Eas.(|75|) and 
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H77|) contain factors of a^~'^ which are 1 or 1/a, respectively. For the Adler function Im[2?p2i] is 
obtained by making the change Ei — > Ei it ivr in the first hne of Ea. ()3U() for Q^ > A^, and in the 
second hne for Q^ < A^. For continuity of the Im part at Q^ = A^ one needs to choose the sign of 
ITT oppositely in the two regions. One then finds for Q^ > A^ 

• ^2 \ ("+i) 



Im[pg(Q2)] 



zbzTT 
1 



a(Q2) 



/ A^ \ 

°° /^2\ 

Yl ^n+l [^0(^ + 1) + Zn+lBi{n + 1)] ( ^J 



n.+l' 



(78) 



Correspondingly, for Q^ < A-^ one finds 



1 



[n 



l2. 



n-l 



a{Q 



EM 

n=2 

oo 



A2 



n-1 



Zn-iAi{n- 1) 



n=2 



A2 



n-1' 



(79) 



Comparing these expressions with Eas. ()64() and (|65|) one then finds that the imaginary part may 
be written directly in terms of the characteristic function ujT){t), 



lm[V^^^{Q^)] 



27rA2 



2^V 



A2 
Q2 



(g' > A^) 



im[pg(g^)] = ±f ^-r (^) (Q^ < A^ 



) 



(80) 



Continuity at Q^ = A^ then follows from continuity of LL'(t) at t = 1. The Cn , and Cn coefficients 
of the OPE and the modified NP expansion are in principle independent of the imaginary part, 
but continuity at Q^ = A^ is dependent upon relations between the Aq^i and So,i residues, such as 
Eas. (|38|) and (|4U|) . and the more complicated structures of Eas. (|69|) - (|7ip . needed for finiteness of 
the Q^ derivatives. Although not strictly necessary for continuity, this continuity follows naturally 
if we write 

K± — 



-D^ikQ^) 



dt uJv{t) + 1 



dujv{t) 



dt 



(81) 



Here k is an undetermined overall real, non-perturbative factor. The t integration here reproduces 
the expressions of Eq.ljHOJ in the two Q^ regions. If the PT component is PV regulated one averages 
over the it possibilities, and combining Ea. H81() with Ea. (|47() for 'Dpj.{Q'^) one can write down a 
result for 'D^^\Q'^) for all values of Q^, 



P(^)(Q2 



dt 



ujv{t)a{tQ^) + K ( ivv{t) + *^^^ ) ^(A' - tQ^) 



(82) 



The Q2 evolution is fixed by the non-perturbative constant k, and by A. The infrared limit is 
^^■^^(O) = 0, we have already noted that 'Dpj,{{)) = 0, the NP component also freezes to zero since 
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FIG. 6: The bold curves corresponding to k = are the perturbative corrections to the observables in 
Eas. H2U|l - H23|l . as in Fig. ^ The upper and lower curves correspond to the overall result including NP 
contributions with the choice k — 1 and k — —1, respectively. 

on integrating the second term one finds an IR limit of w^ (1) — uj^ {1) = 0, from continuity of 
the characteristic function at t = 1. The same expression holds for the other Euclidean observables 
/Cp^((5^) and Upj,{Q'^) on replacing uJt){t) by uJK.{t) and u>u{t), respectively. We plot in Fig. [HI 
the overall result for V^^\Q^), KL^^^Q"^) and U^^\Q'^) for the choices k = 0, k = 1 and «; = -1. 
For the DIS sum rules uJic{t), ^u{t) and their first derivatives are continuous at t = 1. In the 
case of lAj^p{Q'^) there are a total of three non-perturbative terms, and hence the continuity of 
the characteristic function and its first derivative fixes the form of the function up to an overall 
constant factor. Thus Eq. ()82|) does indeed hold for IA^^'{Q'^) without conjecturing the form of 
Eq. (jHU). 



VI. 



INFRARED FREEZING BEHAVIOUR OF R^^-^- 



We turn in this section to a consideration of freezing behaviour of the Minkowskian quantity 

n 

Re+e- which was discussed in Ref.^]. This treatment was criticised in Ref.^], which argued that 
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in fact there is an unphysical divergence in the infrared hmit. We wish to address these criticisms. 
Rf.+e-{s) will be defined by Ea. (|2UI) with the perturbative corrections V{Q^) replaced by 7^(s), -y/i 
here is the e^e~ cm. energy. lZ{s) is related to T){—s) by analytical continuation from Euclidean 
to Minkowskian. One may write the dispersion relation 

If V{t) is represented by a Borel representation as in Eq.@ one arrives at 

Kg(,) = yf~<i«-/«C.5HMB[z,g](., ^ (84) 

There is now an extra oscillatory factor of sin(7r6z/2)/(7r6z/2) arising from the analytical continu- 
ation. In consequence each individual IR or UV renormalon contribution at Q^ = A^ will be finite, 
and the cancellation of Ea. p9() is not required. One can also analytically continue the one-chain 
skeleton expansion result for 'Dprf{Q'^) to obtain 

7^g(s) = 1- Tdt ujv{t) arctan M^^\ . (85) 

Here the principal branch of arctan is assumed so it lies in the interval [— 7r/2, -|-7r/2], and 
arctan(O) = 0. This form is equivalent to the Borel representation of Eq. (|84|) for s > A^, and 
to the modified Borel representation for s < A^. Notice that the choice of principal branch is 
crucial if the PV Borel sum is to be continuous at s = A'^ . The result freezes to the IR limit 
TZpf,[0) = 0, since arctan(O) = on the principal branch. This freezing limit differs from that 
found in the APT approach [9|, where a freezing to an IR limit of 2/6 occurs. This freezing limit 
was also erroneously claimed in Ref. 0, but then the PV Borel sum is discontinuous. In Ref.[3] 
unphysical singularities in the region — A^ < s < lead to extra terms and they find 

7^^^(s) = —r dt ujv{t) arctan ( ^ ") ^^ ] + -r dt uJv{t) + - j dt uj^{t) . (86) 

TTO Jo \ ^ J '^ Jo " J-A'^/s 

These extra terms may be treated as contributions to TZ^piQ"^). The final term leads to an infrared 
divergence as s — > 0, and has an expansion of the same form as the OPE. Notice, however, that the 
Minkowskian OPE for TZ{Q^) is pathological and contains delta-functions 6{s) and their derivatives 
25| . It is only when a smearing procedure in Q^ is used j2fl] that it makes sense. In contrast for 
Euclidean quantities the regular OPE is potentially well-defined, and no smearing is required. 

We will now consider the evaluation of the PV Borel integral for TZpj,, and correct the erroneous 
statements made in Ref. J2] noted above. This can be expressed in terms of generalized exponential 
integral functions Ei(n,it;), defined for Re w > hy 

/oo p—wt 
dt^^ ■ (87) 
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One also needs the integral 



dze-^/ 



^sH^hz/2) _^^^^^^(._ha 



(88) 



Here the principal branch of the arctan is again assumed. Care needs to be taken when Re w <Q. 
With the standard continuation one arrives at a function analytic everywhere in the cut complex 
■w-plane, except at u; = , and with a branch cut running along the negative real semi-axis. 
Exphcitly |22| 



Ei(n,ty) 



i-w) 



n-l 



(n-1)! 



-\n.w 



n—l ^ 

7B+ y] — 

'^ m 



m=l 



m=0 



i-w) 



{m — n+ l)m\ 



(89) 



The \nw contributes the branch cut along the negative real semi- axis. To obtain the PV of the 
Borel integral one needs to compensate for the discontinuity across the branch cut and make 
the replacement Ei(n,t(;) — > Ei(n,t(;) + insignilra w). One then finds that the IR renormalon 
contributions to TZpj^{s) can be written in terms of the functions 

e-^i/''('){-l)nui[Ei{p, -Zg/a{s) - mhzg/2)\ 

-7rRe[((zg/a(s)) + mh/2Y~^] e{s - A^) . 



4>~iP,Q) 



e-^.M^)(_l)9^P-i 



{p-l)l 
The UV renormalon contributions can be written in terms of the functions 



(90) 



cP+{p,q) = e'''/^^'^\-l)nm[Ei{p,zJais)+nTbzj2)] 



+ 






g~p-i 



p-i-i 



-7rRe[{{- zg/a{s)) - iirb/2f-'] e{K' - s) 



(91) 



The PV regulated TZpj^{s) is then given for all values of s by 



n^pks) 



nfT{s)\uv + nf^{s)\iR 



^-arctan 

TTO 



TTba{s] 



\ 2 °° 

+ ^ E iMJ)MhJ) + (Aoij) - A^{j)z, 



+ ^^0-(i,2)+^f;(i?o(j>4i,j) + (5o(j)+i?i(j>, 

i=3 



J' 



(2,j)) ■ (92) 



Note that the presence of the ^-functions is crucial in Eas. (|90() and ()91|). The terms they multiply 
are the extra contributions necessary to obtain the PV when Re w < 0. For s > A^ the second 
contribution is required for the IR renormalon contribution, but for s < A^ it must be switched off, 
otherwise the Borel integral will not be correctly evaluated. With a{s) < for s < A^, Re w < 
occurs for the UV renormalon contributions and the extra term must be switched on to obtain a PV 
regulation of the UV component. Leaving out the ^-function in Eq.(90) would cause an unphysical 
divergence in the infrared, and leaving it out in Eq.(91) would cause asymptotic freedom to fail in 
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the ultraviolet. If the PV is correctly evaluated with arctan remaining on the principal branch for 
s^ < A^ then one obtains 7lpj,{0) = (2/7r6) arctan(O) = 0. Notice that at first sight the PV result 
appears to be discontinuous at s = A^, as the ^-function contributions switch over. However the 
discontinuity is given by the 4'±{l,j) terms, and one finds, upon summing them, a discontinuity 

I f2[Bo{j){-iy + AoU){-iy] = ^i?o(2) = I . (93) 

Here the relation of Ea. (|4U() ensures pairwise cancellations of terms, and Bq{2) = 1 is left over. 
If we remain on the principal branch, however, the arctan term has an equal discontinuity of 
opposite sign, since (2/7r6) arctan (+(X)) = 1/b, whereas (2/7r6) arctan(— cxd) = —1/6, and overall 
there is continuity at s = A^. In Ref.[4] it was wrongly claimed that the PV result is discontinuous 
at s = A^, and instead it was suggested to use a regulation where one throws away the second 
terms in Eqs.(90) and (91). These terms are of the form (A^/s) , and (s/A^) , respectively, and 
so they can simply be absorbed into the regular OPE and its modified form. 

We finally discuss the ambiguous Im[7^p^(s)]. This may be straightforwardly evaluated as 
Im[7^i,^^(s)] = ±ivr^[i?o(n) + i?i(n)z„]z„(-l)"(^J {Q^ > A^) 

oo / A 2 \ '^ 

lm[n%{s)] = TinY^iMn) - A^in)znWi-ir( — ] {Q^ < A^) . (94) 

If one defines u)x>{t) = u)jy {t)+\n.t w^ (t), the split being into the single and double pole renormalon 
contributions, then comparing with Eqs. (|(i4j) and (|fi5j) one finds 

Im[KS(.)l = ±^^c.g"« (^) (. > A=) 

Im[Rg(.)l = ±~^T" (=^) (' < A') ■ (95) 

Notice that only the double poles contribute since the sm{TTbz /2) / (irbz /2) analytical continuation 
term in Eq. (jHSj) contains zeros at z = ibz„ which nullify the single pole contributions. Whilst 

(2) 

the characteristic function ijJx>{t) is continuous at t = 1, the uj)-, {t) function is discontinuous at 
t = —1. The discontinuity is ±2/6 and arises from the same sum in Ea. H93|) which gives an apparent 
discontinuity in the PV TZpj,{s) component, although in the PT case this is cancelled by the arctan 
term. Thus defined in this way lin.[TZ pj.{s)] is discontinuous at s = A^. It would seem that the 
proper way to proceed is rather to use the dispersion relation of Ea. H83|) to analytically continue 
into the Minkowskian region the expression for V^^'Q^ arrived at in Eo. 1)82(1 . Unfortunately the 
one-chain skeleton expansion form for 'D{Q'^) is hard to consistently analytically continue, which 
was a key motivation for the alternative inverse Mellin representation introduced in Eef . J28l| . We 
shall defer further discussion of the more subtle issue of Minkowskian freezing until a later work. 
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VII. DISCUSSION AND CONCLUSIONS 

We have shown in this paper that in the approximation of the one-chain QCD skeleton 
expansion (leading-6 approximation), the perturbative corrections to the parton model result 
for Euclidean observables undergo a smooth freezing to an infrared limit of zero. We explicitly 
studied the Adler function, GLS sum rule and polarised and unpolarised Bjorken DIS sum rules as 
explicit examples, and found that they changed sign in the vicinity of Q^ = t^, and then froze to 
zero at Q^ = 0. Continuity and finiteness at Q^ = Is? follow from continuity of the characteristic 
function u){t)^ and its derivatives at t = 1. The one-chain term is equivalent to the standard 
Borel representation of Eq.(|Sl for Q^ > A'^, and to the modified Borel representation of Ea. H12() . 
previously proposed in Ref.|2|, for Q^ < A^. For the Adler function we established a dictionary 
between the residues of the IR and UV renormalon singularities, and the series expansion 
coefficients of uj{t). Continuity of i^v{t) and its first three derivatives at t = 1 implies relations 
between the residues of IR and UV renormalon singularities: Eqs. (|HH|) . (|in|) . (|H!I|) . (|7n|) and (|7T|) . 
IR renormalons for Q^ > I? lie on the contour of integration in the Borel representation, and 
similarly UV renormalons lie on the contour of integration in the modified Borel representation for 
q2 ^ ^2_ jj^ both cases these singularities lead to an ambiguous imaginary part in T)pr^{Q'^), which 
can be cancelled against an ambiguous imaginary part in the coefficients of the non-perturbative 
terms, in the two Q^ regions. The ambiguous imaginary part may be written directly in terms 
of the characteristic function, as in Ea.H80|). and is continuous at Q^ = A^. If the real parts of 
the condensates are to result in a T)^p{Q'^) which is continuous at Q^ = A^ this suggests that 
one should write these in terms of the characteristic function as well, which led us to conjecture 
Ea. (|81() in which there is a real overall non-perturbative factor k which is undetermined and 
observable-dependent. All of these properties and results hold in general for Euclidean observables 
for which a one-chain result of the form of Eq.Q can be written down. As pointed out in Ref.|j| 
this is not possible for Minkowskian observables such as R^+e- > and in this case the question of 
freezing is more delicate. There is no characteristic function uj'ji{t) to underwrite the smooth 
transition through s = A^ from UV to IR. Indeed the Minkowskian gluon condensate OPE 



contribution is proportional to 5 (s) 23] , so without a smearing procedure it will give an apparent 



infrared divergence as s — > 0. We corrected some erroneous statements about the continuity of 
TZprj,{s) at s = A^ made in |2]. The issue of Minkowskian freezing is interesting and requires 
further investigation. 

An interesting feature of the skeleton expansion representation of Eq. (jH)) concerns the definition 
of a(Q^) for Q^ < A"^. In this region the result of Eq.@ is not in fact the solution of the 
RG equation, but is an analytical continuation of the Q^ > A^ result. However notice that 
in evaluating 'Dpj.{Q'^) for Q^ > h? one is integrating over the region t < I? /Q"^ where the 
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analytically continued a{Q'^) is required. Also notice that QED skeleton expansion results can be 
obtained simply by interchanging UV and IR renormalons, and the limits Q"^ ^> co and Q^ -^ 0. 
The implication is that QCD in the IR energy region is QED-like, and conversely that QED in the 
UV energy region beyond the Landau ghost is QCD-like. Our result of Eq.((S2) is closely related 
to successful models for power corrections based on isolating the IR renormalon ambiguity, such as 
[29|, and to the power correction model of |30l.l3ll|. The latter postulates an infrared finite running 
coupling, and uses a dispersive approach. The infrared limit of the coupling ZSq is a universal 
parameter in this picture, whereas k in our approach is expected to be observable-dependent. A 



more sophisticated discussion of power corrections to DIS sum rules has recently appeared in 32 1 . 
Notice that continuity at Q^ = A^ is the key constraint leading us to suggest Eq. (|81|) . which 
is arguably not a model for power corrections but the actual form of the NP component in the 
one-chain (leading-6) approximation. In future work we intend to report on fits of k and A to 
experimental data on the Adler function (e.g., the analysis of |2d|)) and DIS sum rules. It would 
also be interesting to compare our results for the Q^-dependence of the polarised and unpolarised 
Bjorken sum rules in the light of the relations between them noted in Ref.|3J]. Whilst the freezing 
is straightforward to analyse in the leading-6, one-chain approximation, it is much harder to analyse 
at the two-chain level, where the anomalous dimensions for the operators will enter. The IR^^UV 
conformal relations between renormalon residues and condensate coefficients would appear to be 
the key ingredient in the freezing picture, and there is hope that they can continue to hold at higher 
orders in the QCD skeleton expansion, if indeed such an expansion can be consistently formulated 



^ 



There is clearly much still to investigate. 
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